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Vh ' Abstract 



We argue that the "vortex-finding" property of maximal center gauge, i.e. the ability 
of this gauge to locate center vortices inserted by hand on any given lattice, is the key 
to its success in extracting the vortex content of thermalized lattice configurations. We 
explain how this property comes about, and why it is expected not only in maximal center 
gauge, but also in an infinite class of gauge conditions based on adjoint-representation 
link variables. In principle, the vortex-finding property can be foiled by Gribov copies. 
This fact is relevant to a gauge-fixing procedure devised by Kovacs and Tomboulis, where 
we show that the loss of center dominance, found in their procedure, is explained by a 
corresponding loss of the vortex-finding property. The association of center dominance 
with the vortex-finding property is demonstrated numerically in a number of other gauges. 



1 Introduction 

Numerical evidence in favor of the center vortex theory of confinement has been steadily 
accumulating over the past three years ||l|-0] . Underlying most of these numerical studies 
is a technique for locating center vortices in thermalized lattice configurations, known as 
center projection in maximal center gauge. 

In its "direct" version |@,^, maximal center gauge is the gauge in which 

R = '^'^\T^^[U^{x)] is a maximum (1) 

X fj, 

This gauge brings each link variable as close as possible, on average, to a Zn center element, 
while preserving a residual Zn gauge invariance. Center projection is a mapping of each 
SU(A^) link variable to the closest Z^ center element; e.g. in SU(2) gauge theory, center 
projection is the mapping 

U^ix) -^ Z^{x) = signTr[f/^(a;)] (2) 

The excitations on the projected Zn lattice are point-like, line-like, or surface-like objects, 
in D = 2, 3, or 4 dimensions respectively, known as "P-vortices." These are thin objects, 
only one lattice spacing across. There is substantial numerical evidence, for SU(2) gauge 
theory, that thin P-vortices lie roughly in the middle of thick center vortices on the un- 
projected lattice, and that these thick vortices produce the entire asymptotic SU(2) string 
tension |^ . The number of P-vortices mod 2 linking a large loop is closely correlated with 
the sign of the corresponding SU(2) Wilson loop [0,0, and P-vortices themselves are regions 
of high action on the unprojected lattice 0. It is found that removal of center vortices 
not only removes the asymptotic string tension, but chiral symmetry goes as well, and the 
SU(2) lattice is then brought to trivial topology 0. The vortex density has been found 
to scale as predicted by asymptotic freedom 0, 0,|[, and, at finite temperature, the non- 
vanishing string tension of spatial Wilson loops in the deconfined phase can be understood 
in terms of vortices winding through the periodic time direction [0,P|. Vortex percolation 



properties at finite temperature have also been studied in ref. |T0|. The world-lines of 



abelian-projection monopoles are found to lie on P-vortex surfaces, and the field-strength 



associated with these monopoles seems to be collimated in the vortex direction ||T2[. Fi- 
nally, it appears that even the Casimir scaling of higher-representation string-tensions at 
intermediate distance scales can be understood in terms of the finite thickness of center 
vortices p. 

On these grounds, we are confident that the vortices identified by our gauge-fixing + 
projection procedure are physical objects which are crucial to the confinement mechanism. 
But a disquieting question remains, namely: Why does this procedure work? In what way 
does the gauge choice (|I|), combined with the projection (||), identify center vortices? In 
particular, since a vortex creation operator (unlike a monopole creation operator) makes 
no reference whatever to any special gauge choice or Higgs field, why do we need to fix to a 
definite gauge in order to locate vortices? These questions become quite urgent when it is 



recognized that apparently minor changes in the gauge-fixing condition, or even, as shown 



recently by Kovacs and Tomboulis |]13|, a small change in the gauge- fixing procedure^ can 
be catastrophic, and the resulting P-vortices no longer correspond to anything physical. 
So what crucial property of the gauge-fixing/projection procedure has been lost, when the 
method fails? 

In this article we identify this crucial property of our procedure as the "vortex-finding" 
(VF) property, by which we mean the following: Suppose, in any given thermalized lattice, 
a center vortex is inserted "by hand" via a discontinuous gauge transformation. The lattice 
now contains at least one center vortex in a known location. Upon gauge-fixing and center 
projection, a set of P-vortices is identified. Is the vortex inserted by hand found among this 
set of P-vortices? If so, then the procedure has the vortex-finding property. This property 
seems like a reasonable demand to make of any method which is advertised to extract the 
vortex content of lattice configurations, and is presumably a necessary condition for its 
success. 



2 The VF-Property in Adjoint Gauges 

Let us examine why maximal center gauge and, as we will see, an infinite class of other 
gauges, might have this vortex-finding property. We begin by noting that the maximal 
center gauge defined by eq. (0) is really the adjoint Landau gauge; i.e. it is equivalent to a 
Landau gauge-fixing condition on adjoint links 

i? = ^ ^ Tt[Ua^i{x)] is a maximum (3) 

X n 

where Uaii{x) is the link variable in the adjoint representation. This motivates considering 
other gauge-fixing conditions of the general form 

7^[t/A] is a maximum (4) 

such that the gauge condition 

1. depends only on the adjoint representation links; 

2. is a complete gauge-fixing of the adjoint link variables; 

3. transforms most links to be close to center elements, at weak coupling. 

We will refer to these as "adjoint" gauges. 

Let U denote some thermalized lattice configuration. A center vortex is created, on 
the background U, by a discontinuous gauge transformation (some explicit examples will 
be given below). Denote the resulting configuration as U', whose field strength differs 
from that of U only at the vortex core. Away from the vortex core, the corresponding 
link variables of each configuration in the adjoint representation, denoted Ua and U'a, are 



gauge-equivalent. This is because the global discontinuity of the vortex-creating gauge 
transformation, given by a center element, is invisible in the adjoint representation. 

The crux of the argument is this: It is assumed that the gauge-fixing condition @) is 
complete for adjoint links. If we for the moment ignore both (i) the Gribov copy problem; 
and (ii) the region of the lattice corresponding to the core of the created vortex, then Ua 
and [/^ are gauge equivalent, and gauge-fixing to (^) should map both Ua and U'j^ into 
the same gauge- fixed configuration Ua- Under these mappings, the link variables U and 
U' in the fundamental representation are transformed to configurations U and U' corre- 
sponding to the same adjoint configuration Ua- Since they correspond to the same S0(3) 
configuration, the U and U' lattice SU(2) configurations can differ only by continuous 
and/or discontinuous Z2 transformations. Then, because a continuous gauge transforma- 
tion, associated with the gauge-fixing, cannot undo the discontinuous transformation which 
created the center vortex, the vortex originally inserted in U' appears as a discontinuous 
Z2 transformation relating U' to If. 

What has happened here is that the original discontinuous gauge transformation, which 
may be quite smooth (up to the discontinuity) and extended, has been squeezed by the 
gauge-fixing condition to the identity everywhere except on a Dirac volume (bounded by 
the vortex core), where it has the effect of simply multiplying a certain set of links by —1. 
Upon center projection, U —>■ Z and U' -^ Z' , and the projected configurations differ by 
the same discontinuous Z2 gauge transformation. This discontinuity then shows up as an 
additional P-vortex in Z', not present in Z, at the location of the vortex inserted by hand. 

This "vortex-finding property" goes a long way towards explaining the success of maxi- 
mal center gauge in extracting the vortex content of lattice configurations. If, in fact, lattice 
vacuum configurations have the form of a product of vortex-creation operators which op- 
erate on a non-confining background state, then a procedure with the VF-property may 
be reasonably expected to locate these confining vortices. The above argument not only 
explains why maximal center gauge should have the VF-property, but also suggests that 
any complete gauge-fixing condition on adjoint links might have the VF-property as well. 
However, there are two ways that the argument we have presented may go wrong: 

Gribov Copies: The argument for vortex-finding is based on complete adjoint gauge- 
fixing; i.e. any two gauge-equivalent SO (3) configurations should be mapped to a 
unique adjoint link configuration. Unfortunately, standard methods of implementing 
eq. (^ in practice, such as over-relaxation and simulated annealing, usually wind up 
in local maxima of 7?.[?7a], known as Gribov copies, rather than the global maximum. 
So the argument for the VF-property can fail at this point. 

Vortex Cores: The configurations Ua and t/^ are only gauge-equivalent outside the 
vortex cores. Because they are gauge-equivalent outside this (relatively) small region, 
we have assumed that Ua and U'j^ will transform to the same gauge-fixed SO (3) 
configuration Ua outside the core region. This assumption, however, could simply be 
wrong. 



Because of these caveats, we have no proof of the vortex-finding property in the adjoint 
gauges. This is just as well, since we will soon discuss cases where the property fails. 
Nevertheless, we now have some idea of why the VF-property might hold in maximal 
center gauge, and a motivation to test this property numerically, both in maximal center 
and in other adjoint gauges, to see if it is correlated with center dominance. We should 
also note that our argument for the VF-property has no loopholes when (i) the inserted 
vortex is thin (one lattice spacing thick), so that Ua and f/^ are SO (3) gauge-equivalent 
everywhere on the lattice; and (ii) there is no Gribov problem. In that case, according to 
our argument, center projection is certain to find the inserted vortex, as will be confirmed 
in subsection 3.8 below. 

3 Testing the VF-Property 

We begin by describing a class of discontinuous SU(2) gauge transformations gv on an 
L^ xT lattice, which create two parallel thin vortex surfaces at time t = T that are closed 
by lattice periodicity. These transformations have the usual form 

U^{x) -> U'Jx) = gv{x)Uf,{x)gl,{x + jl) (5) 



except that the gauge transformation has the discontinuity 



lias LliC UlSUUliLlliUlLJ' 

{ -gv{x,y,z,l) xi<x<X2 
y gv{x,y, z^l) otherwise 



Its not hard to see that each discontinuous transformation gy of this form is equivalent to 
the mapping 

U4,{x,y,z,T) ^ -U4,{x,y,z,T) for xe[xi,X2] (7) 

with all other links unchanged, followed by an ordinary continuous gauge transformation 

M^) ^ K(^) = g{x)U^{x)g\x + fi) (8) 

Suppose that the transformation by gy is performed on a thermalized lattice, and 
that the configurations U and U' are gauge-fixed to an adjoint gauge, and then center- 
projected. The corresponding projected configurations are denoted Z and Z', and we 
denote the Polyakov lines in these projected configurations by P{x, y, z) and P'(x, ?/, z) 
respectively. If the gauge-fixing + projection procedure has the vortex-finding property, 
then the Z' lattice should contain P-vortex surfaces on the dual lattice at fixed t = T and 
fixed X = Xi — 1,X2. The two parallel surfaces bound a Dirac 3-volume, e.g. as indicated 
in fig. |1] (apart from its boundary, the location of the Dirac volume is gauge- dependent). 
Now, if all the other P- vortices were located in identical positions in Z and Z' (i.e. if, apart 
from the inserted vortex, Z and Z' were equivalent up to a Z2 gauge transformation), then 
the vortex-finding property is verified if 
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Figure 1: Left: The effect of the transformation (J^) on hnks at fixed y,z and t = T, 1. 
Downward pointing time-hke hnks have been multiphed by —1. Right: Parallel vortex 
surfaces created at constant xi,X2,t = T bound a Dirac volume, closed in the y, z directions 
by lattice periodicity. 

At this point, however, we have to face up to the Gribov copy problem. Center vortices 
on unprojected lattice configurations are comparatively thick (f» 1 fm) objects, and the 
precise "middle" of a vortex core, which P-vortices are supposed to locate, is a little 
ambiguous. In fact, it is found that if we take two gauge copies Ui and Un of the same 
configuration ?7, gauge-fix each to maximal center gauge and then center project, that the 
P-vortices in the corresponding projected configurations Zj and Z//, although correlated, 
differ somewhat in position. Thus, a P-vortex plaquette in Zi which is just inside the 
minimal surface of a given loop C, may lie just outside this minimal surface in configuration 
Zii (for a more detailed discussion and related numerical results, see ref. P|). 

The randomizing effect of small differences in P-vortex locations in Z and Z' (which cor- 
respond to the same "thick" vortices in f/ and f/') will cause the vev {P'{x, y, z)P{x, y, z)) 
to differ from —1 and -|-1, respectively, inside and outside the region x G [a;i,a;2]. However, 
if the vortex-finding property is valid, then eq. (^ should hold true when this random- 
izing effect is factored out. To this end, we generate a set of thermalized SU(2) lattice 
configurations and, from each configuration, we obtain three configurations consisting of 

I. The original configuration, denoted Ui. 

II. The original configuration with an inserted vortex, denoted f///. In most of our 
simulations, a thin vortex is inserted by applying the mapping (|^), followed by a 
random (but continuous) gauge transformation gu at every site. 

III. A gauge copy of f//, denoted Uui, obtained by applying a random continuous gauge 

transformation gm to Uj at every site. 



We then fix each of these configurations to an appropriate adjoint gauge, center project 



according eq. (|]), and compute the ratio 

^ ' j:y,APi{x,y,z)Pjjjix,y,z)) ^ > 

where P/, P//, Pm denote Polyakov hnes in the corresponding projected configurations, 
and where the denominator has the effect of factoring out the randomizing (Gribov-copy) 
effects just mentioned. Then the gauge-fixing + projection procedure has the vortex-finding 
property only if 

1+1 otherwise 

3.1 Thin Vortex Insertion, Maximal Center Gauge 

For our first test of the VF-property, we have performed the Monte Carlo simulation on 
a 14^ X 12 lattice at /5 = 2.3, and applied a discontinuous gauge transformation (eq. (0) 
followed by a random continuous gauge transformation) to insert thin vortices as described 
above. We then fix to maximal center gauge via over-relaxation, as described in ref. [Q, 
center project according to eq. (H), and calculate G{x). The discontinuity in eq. (|]) is 
chosen to lie in the range x G [4, 10]. 

The result for G{x) is shown in fig. |^. The criterion for the vortex-finding property, eq. 
(|rT|), appears to be nicely confirmed in this case. Given the successes of our procedure, 
summarized in the Introduction, in extracting the vortex content of vacuum configurations, 
the existence of a vortex-finding property was perhaps to be expected. We find that in 
this case, the existence of Gribov copies in our gauge-fixing procedure does not destroy the 
VF-property. 

Instead of inserting vortices via the transformation (J^), followed by a random gauge 
transformation (H), we have also considered vortex insertion by gauge transformations 
which are smooth in the t-direction up to the discontinuity at t = T + 1, taking the form 

, . j exp in{t - 1)(73/T for xe[xi,X2] ,-.r,^ 

gv{x,y,z,t) = i L ^ ^ ' J _ l (•^2) 

[ 1 otherwise 

and performed simulations for the same parameters as before (/3 = 2.3, 14^ x 12 lattice, 
[xi, X2] = [4, 10]). As in the previous case, the results are consistent with the vortex- finding 
property in eq. (0). 

Apart from our tests here, based on the values of G{x), we would also like to mention 
some relevant results reported recently by Montero in ref. [H]. Montero, building on the 



work of ref. |T5[, constructs classical SU(3) center vortex solutions on a periodic lattice. 
The stability of the solution is due to the use of twisted boundary conditions, with two 
dimensions of the periodic lattice chosen much smaller than the two other dimensions. The 
lattice is fixed to maximal center gauge and then center-projected. It is found that P-vortex 
plaquettes accurately locate the middle of the classical vortex. This is very interesting 
independent evidence for the existence of the vortex-finding property of maximal center 
gauge (in SU(3) gauge theory, in this case). 



Direct to MCG: T=12, beta=2.3 
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Figure 2: Graph of G{x) for configurations with thin inserted vortices (14^ x 12 lattice, 
380 lattices, (3 = 2.3). The gauge discontinuity (|^) is located in the 3- volume at x G [4, 10], 
t = T. 



3.2 Thin Vortex Insertion, Kovacs-Tomboulis Procedure 

The argument for the vortex-finding property of adjoint gauges, as presented in section 
2, explicitly neglects the problem of Gribov copies, and in fact the possibility of gauge- 
fixing to a unique configuration was an essential step in the argument. Conversely, it 
follows that the existence of Gribov copies in the gauge fixing procedure might destroy 
the vortex-finding property, and, as a consequence, the ability of the projection procedure 
to locate center vortices in thermalized lattice configurations. This appears to be what 
happens in a modification of our gauge-fixing procedure considered recently by Kovacs 
and Tomboulis [jl3| (and also noted in ref. 0.) 

Kovacs and Tomboulis suggest fixing first to the (usual) lattice Landau gauge, before 
fixing to maximal center gauge via over-relaxation. Of course, if the over-relaxation proce- 
dure transformed each configuration to a unique, global maximum of eq. (|ip (i.e. if there 
were no Gribov copies), then a preliminary fixing to lattice Landau gauge could not pos- 
sibly make a difference to the result. However, the over- relaxation procedure only finds a 
local maximum of (|I|). Thus the starting point can make a difference, at least in principle.^] 
Kovacs and Tomboulis find that with a Landau gauge starting point, center dominance is 
completely lost in the projected configurations, which appear to have no asymptotic string 



^There is not much difference, fiowever, in R. Making a set of gauge copies of a given configuration, 
and gauge-fixing each copy to maximal center gauge either directly, or with a Landau-gauge starting point, 
we find that the variation in R of eq. (0) among different copies far exceeds the average difference in R 
between the two procedures; in fact, our own results for this average difference are not yet statistically 
significant. 



tension. This fact does not really affect the status of P-vortices, identified by the usual 
method, as locators of physical objects. That status is well established by the correla- 
tion of P-vortices with gauge-invariant observables, and by their scaling properties. In the 
modified procedure, however, it would appear that some important property, essential for 
extracting the center vortex content of vacuum configurations, has been lost. From every- 
thing that has been said until now, the obvious candidate for this "essential property" is 
the vortex-finding property. 



Landau Before MCG: T=12, beta=2.3 
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Figure 3: Graph of G{x) for configurations with a thin inserted vortex. Configurations are 
first fixed to the Landau gauge, and only then to maximal center gauge (14^ x 12 lattice, 
260 lattices, (3 = 2.3). The Dirac volume is the same as in fig. ^. 



In fig. ^ we display our results for the observable G{x), in which each configuration 
Uj, Uji, Ujji is first fixed to Landau gauge, before fixing to maximal center gauge via over- 
relaxation. The contrast between Figs. ^ and |^ is quite striking; the Kovacs-Tomboulis 
procedure is clearly inconsistent with eq. ([TTl) , and seems to have completely lost the vortex- 
finding property. Only at the boundary of the Dirac volume, where there is a strong local 
field strength, does G{x) show some effect (although G{x) > even there). But in the 
middle of the Dirac volume, G{x) is comparable to its value outside the volume. It seems 
likely that if the interval [xi, 0:2] were large enough, the middle of the Dirac volume would 
be indistinguishable from the outside region. Since this modified procedure cannot even 
identify a thin center vortex inserted into the lattice by hand, there is no reason to expect 
it to locate the fuzzy, much more diffuse center vortices generated by the gauge theory 
dynamics. 



3.3 Thin Vortex Insertion, Asymmetric Adjoint Gauge 

The argument for the vortex-finding property in section 2 did not single out the maximal 
center gauge in particular; it is possible that other adjoint gauge choices might work just 
as well. Let us therefore introduce the "Asymmetric Adjoint Gauge" 



R = Y.Y.^Mu,ix)] 



IS a maximum 



(13) 



X fl 



where c^ is some set of four positive numbers. For an exploratory run aX f3 = 2.3, again on 
a 14^ X 12 lattice, we chose a set of values 

{ci, C2, C3, c4 = {1.0, 1.5, 0.75, 1.0} (14) 

and fixed to the gauge (|1|) by over- relaxation (without prior Landau gauge- fixing) . 
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Asymmetric Adjoint Gauge: T=12, beta=2.3 
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Figure 4: Graph of G{x) for configurations with a thin inserted vortex. Configurations 
are fixed to the asymmetric adjoint gauge with c^ = {1, 1.5,0.75, 1} (14^ x 12 lattice, 380 
lattices, f3 = 2.3). The Dirac volume is the same as in fig. Q. 

The result for G{x) in this new gauge is shown in fig. ^. Once again, the vortex-finding 
property is satisfied within errorbars. The next question is: Do the corresponding center- 
projected configurations display center dominance (i.e. do they have the same string tension 
as the unprojected configuration)? In a run at /5 = 2.3 on a 14^ lattice, in the asymmetric 
adjoint gauge ([T3|) , ([Mf) , the projected Creutz ratios do, in fact, agree quite well with 
the asymptotic string tension extracted from gauge-invariant Wilson loops, reported in 
ref. [0. The results are shown in fig. §. 
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Center Projected Creutz ratios, 14 lattice, p=2.3 
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Figure 5: Center dominance for the gauge with c^ = {1, 1.5, 0.75, 1}. (14^ lattice, /3 
Sohd hne shows the value of asymptotic string tension quoted in ref. [|T^ . 



2.3.) 



3.4 The Modulus Landau and Adjoint Coulomb Gauges 

It was certainly not obvious that the preliminary Landau gauge-fixing, suggested by Kovacs 
and Tomboulis, followed by maximal center gauge-fixing via over-relaxation, would destroy 
the VF-property (and, as a consequence, center dominance) in the projected configurations. 
It is also a little surprising that certain choices of adjoint gauge, which would naively seem 
just as good as maximal center gauge, appear to have similar problems. An example is a 
"modulus" version of the usual lattice Landau gauge 



^ = EE|Tr[t/,(x)] 



IS a maximum 



(15) 



which is also an adjoint gauge, as defined by the conditions set out in section 2. We 
believe, on the basis of the argument in section 2, that if R could be fixed to a unique 
global maximum then this gauge would also have the VF-property, and the projected 
configurations would exhibit center dominance. However, as with maximal center gauge, 
the only known gauge-fixing techniques are simulated annealing and over-relaxation, and 
these are both plagued with Gribov copies. So the only way of testing the vortex-finding 
property, and center dominance, is numerically. 

We have used over-relaxation at /3 = 2.3 to fix to the modulus Landau gauge. In this 
case, the falloff with area of the projected Wilson loops appears to be much faster than 
that of the unprojected loops, as seen in fig. ^; center dominance is clearly lost. 

The corresponding values for G{x) in modulus Landau gauge have very large errorbars, 
and this is simply because both the numerator and denominator in eq. ([T0| ) have values 
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Figure 6: Center-projected Wilson loops vs. Area in modulus Landau gauge, for square 
Rx R and rectangular Rx (R + l) loops, at /3 = 2.3 on a 12^ lattice. Unprojected Wilson 
loop values are also shown for comparison. 



Modulus Landau Gauge: T=12, beta=2.3 
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Figure 7: Numerator N{x) (eq. (|T6D) of the ratio G{x), for modulus Landau gauge (14^ x 12 
lattice, 380 lattices, /3 = 2.3). The Dirac volume is the same as in fig. ^j. 



which are, within errorbars, consistent with zero. Results for the product in the numerator 






(16) 



12 



are shown in fig. |^, where L (= 14) is the lattice length in the y, 2;-directions. It is clear 
that there is in this case a total loss of the vortex-finding property. In contrast to the 
Kovacs-Tomboulis case, in which G{x) was positive irrespective of linking to the inserted 
vortex, in this case G{x) seems to be essentially 0/0, irrespective of linking number. But 
just as in the Kovacs-Tomboulis case, the failure of center dominance in the projected 
configurations is associated with a corresponding loss of the vortex-finding property. 
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Figure 8: Graph of G{x) for configurations with thin inserted vortices (14^ x 12 lattice, 
110 lattices, [3 = 2.3) in modulus Landau gauge. In this case, the lattice is first fixed to 
maximal center gauge, before fixing via over-relaxation to modulus Landau gauge. 



Since the modulus Landau gauge should be a perfectly good adjoint gauge, the failure 
of the vortex-finding property in this case must be attributed to the Gribov copy problem. 
This impression is strengthened when one compares the rms value of Qq, where 



U„ 



a^I + ia ■ a 



(17) 



in both the modulus Landau and maximal center gauges. We find that Oq™* = 0.76 in 
modulus Landau gauge, vs. a^"^^ = 0.86 in maximal center gauge, ai (3 = 2.3, and that is 
a surprisingly large discrepancy, in view of the fact that both gauges are, in some sense, 
trying to bring links close to center elements. This suggests trying out a variation of the 
Kovacs-Tomboulis procedure: Perhaps modulus Landau gauge-fixing could be improved, 
if we first fix to maximal center gauge, and then fix to modulus Landau gauge via over- 
relaxation. 

It turns out that a preliminary gauge-fixing to maximal center gauge restores both 
center dominance and the vortex-finding property. For center dominance, we find that the 
Creutz ratios, at /3 = 2.3, are clustered near x[R,R] ~ 0.14, quite close to the asymptotic 
string tension obtained on unprojected lattices. The result for G{x) is shown in fig. H; it is 
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similar to what we previously found for maximal center gauge and the asymmetric Landau 
gauge. 

Once again, the presence or absence of the vortex-finding property is associated with 
the presence or absence of center dominance in the projected configurations. 



Adjoint Coulomb Gauge: T=12, beta=2.3 
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Figure 9: Numerator N{x) (eq. (|T6|)) of the ratio G{x), for adjoint Coulomb gauge (14^ x 12 
lattice, 180 lattices, /3 = 2.3). The Dirac volume is the same as in fig. ^j. 

In section 3.3 we studied the asymmetric adjoint gauge, with a moderate variation of c^ 
between 0.75 and 1.5. It is also of interest to look at limiting cases of this gauge, where the 
ratios of some of the c^'s tend to zero or infinity. A particular example we have investigated 
is the adjoint Coulomb gauge, in which 

{C1,C2,C3,C4} = {1.0,1.0,1.0,0.0} (18) 

Center dominance is lost in this gauge also, although the disagreement between full and 
projected string tensions is not quite as bad as in modulus Landau gauge without pre- 
conditioning. Creutz ratios in the projected configurations a.t (3 = 2.3 cluster around 
Xproj{R, R) ~ 0.18, whereas the full asymptotic string tension is around aa^ = 0.135. 

As in modulus Landau gauge, the loss of center dominance in adjoint Coulomb gauge 
is accompanied by a breakdown of the vortex-finding property, as seen from a plot of the 
numerator (|TBp in fig. |^. The qualitative difference between adjoint Coulomb gauge, and 
the asymmetric adjoint gauge studied in section 3.3, may be connected with the fact that 
adjoint Coulomb gauge is not really an adjoint gauge, as defined in section 2. Our third 
criterion for adjoint gauges is that link variables should be brought close to center elements. 
What we find instead, at /? = 2.3, is that while the rms value Oq"^* for spacelike links is 
0.89, the corresponding rms value for timelike links is only 0.50 (and this is identical to 
the rms values of the other three components a™*. Thus, while spacelike links do indeed 
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fluctuate near center elements, the timelike links do not, and the "close-to-center" criterion 
is violated. This criterion for adjoint gauges will be further discussed in section 3.6. 

3.5 Thicker Inserted Vortex, Maximal Center Gauge 

All of the inserted vortices so far are thin vortices; the vortex core has a thickness of one 
lattice spacing. On the other hand, the vortices generated dynamically by the gauge theory 
seem to be rather thick objects, with core widths on the order of one fermi. It is therefore 
of interest to see if an inserted vortex with a somewhat thickened core can also be identified 
on the projected lattice. 

In our previous examples, the —1 discontinuity in the vortex-creating transformation 
gv begins and ends abruptly at x = Xi and x = X2- To create a vortex core several lattice 
units thick in the x-direction, we simply make the transition from the —1 discontinuity to 
-|-1 continuity more gradual. This is done by replacing the mapping in eq. (|^ by 



Ua{x, y, z, T) -^ Ua{x, y, z, T) exp[ia{x)a3,] 



(19) 



where a{x) interpolates smoothly from a{x) = outside the Dirac volume, to a{x) = it 
inside the Dirac volume. 

We make the choice for a{x) shown in Table |I] on an 18'^ x 8 lattice, followed by a 
random gauge transformation, to obtain configurations Uu. 
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Table 1: Thick vortex core: a{x) in eq. (p!9|). 

Maximal center gauge-fixing and center projection are carried out, and G{x) is calculated, 
with results shown in fig. [1^ (obtained at /5 = 2.3). Once again the vortex- finding property 
is quite evident, with G{x) interpolating smoothly from +1 to —1 across the vortex core. 

3.6 The "Close-to-Center" Criterion 

All of the gauges we have considered have the property of bringing most link variables 
close to the ±/2 center variables. This was listed as the third criterion of an adjoint gauge 
in section 2, and we should now explain the rationale behind this criterion. 

The picture which is implied by the numerical studies [1 — 11] is that thermalized SU(2) 
lattice configurations have the form 



U = GvoU. 



NC 



(20) 



where Gy is an operator creating center vortices responsible for confinement, while U^c is 
a non-confining lattice background. The operator Gy is a "smoothed" discontinuous gauge 
transformation; i.e. it is has the form of a discontinuous gauge transformation, away from 
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Figure 10: G{x) computed for thicker vortices on an 18^ x 8 lattice, 380 lattices, P = 2.3. 

vortex cores. ^ In an adjoint gauge which has the vortex-finding property, center projection 
should locate the approximate position of vortices created by Gy- 

There is no guarantee, however, that gauge-fixing -|- center projection will not also 
produce, in addition to the P-vortices associated with Gy, a lot of other extraneous P- 
vortices with no physical content, nor is there any obvious reason that these extraneous 
P-vortices should be independent of the choice of adjoint gauge. This is the motivation for 
introducing, in the definition of adjoint gauges, a "close-to-center" requirement. 

In maximal center gauge (or any other adjoint gauge), not all links are close to center 
elements, even at large /3; some links lying in the middle of center vortices must deviate 
very substantially from ±/2. This deviation from the center elements is necessary, in order 
that a projected plaquette can equal —1, while the corresponding gauge- invariant plaquette 
on the unprojected lattice is close to +1 at weak coupling (cf. [|ll]). The aim of the "close- 
to-center" requirement is that links should deviate substantially from center elements only 
where they must do so, i.e. in the cores of vortices created by Gy', elsewhere the gauge- 
fixing condition should force links to fiuctuate in the vicinity of ±/2. This criterion is, of 
course, most obvious in maximal center gauge. 

Let MC denote maximal center gauge, and AG denote another adjoint gauge which (i) 
has the VF-property, and (ii) has most links close to center elements, except in the middle 
of center vortices. We can then argue that, away from vortex interiors, the center- projected 
configurations Zmc and Zag will be Z2 gauge-equivalent. Let t/lfc ^'ag represent a ther- 
malized configuration U fixed to the MC and AG gauges, respectively. We can write 



ri' 



^MC 



U" 

'^ MC 



where 



U" 

'^ MC 



?>ignIi[U'Mc\U'MC 



^If Gv had this form everywhere, then it would produce thin vortices of very high action; instead, the 
regions of high field strength are smoothed out, and thick vortices are created. 
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U'^a = ZagU'^g where U'^a = signTip'^aP'AG (21) 

where the product and sign trace operations are of course performed hnk by hnk.0 Because 
U'j^Q and Uj^Q are gauge-equivalent, we have 

U'ag = N)°^MC 

= izoZMc)i9oU';jc) (22) 

where z is a Z2 gauge transformation, and TT[g] > 0. Away from vortex interiors, most 
hnks in 1/'^,^^ and U'^q are close to ±J2 at large /3, which implies that g o U'Ijq ~ ±/2 are 
also close to center elements. Combining this with the fact that U'Ijq ~ +I2 and Tr[g] > 0, 
we deduce that g{x) ^ g{x + fi) in this region, and therefore g o U'^^^ ~ +/2. We can then 
identify 

T Til T Til 

Uag = 9°Umc 

Zag = zoZmc (23) 

According to this argument, Zag and Zmc are Z2 gauge-equivalent, away from the vortex 
interior. 

Inside a vortex core, U'Iiq is not necessarily close to the identity, so Zag and Zmc are 
not necessarily gauge-equivalent. This implies that there will in general be some variation 
in P-vortex location among different adjoint gauges (and among Gribov copies in the same 
adjoint gauge), associated with the finite width of the center vortex core. But apart 
from this local variation, the vortices found by two adjoint gauges should be basically the 
same, providing that both gauges have the vortex-finding property, and providing that 
the links in both gauges are close to center elements everywhere except where they must 
strongly deviate, i.e. in the core of vortices created by Gy- How well this "close-to-center" 
requirement is actually fulfilled, in different adjoint gauges with the vortex-finding property, 
is a topic which has not yet been studied in any detail. 

3.7 Speculations about Cooling and Smoothing 

If maximal center gauge + center projection is applied to cooled configurations, then the 
projected string tension is drastically reduced after only a few cooling steps (similar results 
are found for RG-smoothed configurations) 0. But it was also found in ref. |^ that the 
cores of thick center vortices, whose approximate location is identified before cooling, also 
expand considerably after only a few cooling steps, as measured by the one-vortex to 
zero-vortex loop ratio Wi{C)/Wq{C). 

Center vortices in thermalized lattices have been found to be rather thick, "fuzzy" 
objects, which percolate throughout the lattice. It appears that, after a few cooling steps, 
the vortex cores overlap so much that there is virtually no region of the lattice which is 



■^The U'ljfj configuration, incidentally, is the lattice studied by de Forcrand and D'Elia 0, and it is 
found to have neither confining nor chiral-symnietry breaking properties. 

17 



not part of a vortex core. In view of the second caveat ("Vortex Cores") in section 2, it 
is then not so surprising that the maximal center gauge + center projection procedure is 
unable to extract the locations of these very fat, highly overlapping objects (although the 
method seems to have no trouble finding thin, inserted vortices after a few cooling steps). 

3.8 Laplacian Center Gauge 

According to our argument in section 2, an adjoint gauge without the Gribov copy problem 
is guaranteed to have the vortex-finding property, at least for thin inserted vortices where 
vortex width is only one lattice spacing. Such a gauge was invented recently by Alexandrou, 



D'Elia, and de Forcrand |]T^; it is known as the "Laplacian Center Gauge." 

To fix to Laplacian center gauge, one finds the eigenvectors 'u(x), v{x) corresponding to 
the two lowest eigenvalues of of the covariant lattice Laplacian 

K = '^DS^y^"' - E Ut±,i^)5.±,,y (24) 

where 

Ut = InU.crX^'] (25) 

are matrix elements of the link variables in the adjoint representation. The gauge trans- 
formation which takes the lattice configuration into Laplacian center gauge is the trans- 
formation which rotates u{x) to lie in the 3-direction, and v{x) to lie in the 1-3 plane, at 
every site x. This procedure is free of Gribov copies, and it was found in the second refer- 
ence of [|l^ that there is agreement between the asymptotic string tensions on the full and 
projected lattices, at least at the coupling P = 2.4 used there. Although Laplacian center 
gauge does not involve directly maximizing a functional (as in eq. (^)), it does satisfy all 
three conditions below (|), and thus qualifies as an adjoint gauge. 

We have checked the VF-property for Laplacian center gauge as before; i.e. for a thin 
inserted vortex on the same lattice as in fig. ^, at /5 = 2.3. The somewhat startling result 
is that the condition 

G(x) = l -^ ^e[4>10] .26) 

1+1 otherwise 

is satisfied by the numerical data exactly. In other words, there are no errorbars at all; every 
configuration gives the same result for G{x). Fixing to Landau gauge, prior to Laplacian 
center gauge, makes no difference to this result. 

In fact, the absence of errorbars on G{x) in Laplacian center gauge is simply a confirma- 
tion of the reasoning in sections 2 and 3. As noted above, the VF-property is guaranteed in 
this gauge, at least for thin vortices. We have also argued, in section 3, that any deviation 
from eq. (^ must be due to small variations in P-vortex location, from one Gribov copy to 
another, and these variations are responsible for the statistical fluctuations in G{x). Elimi- 
nating Gribov copies eliminates in turn these statistical fluctuations, and the VF-property 
for Laplacian center gauge is seen in the most compelling way. 
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Laplacian center gauge is not quite as "close to center" as maximal center gauge. 
We find, for example, that the rms value of a^ = Tt[U]/2 at /5 = 2.3 is approximately 
^rms _ g gg ^^ maximal center gauge, but only Oq™* = 0.82 in Laplacian center gauge. 



There are also more P- vortex plaquettes in the projected Laplacian configurations; ref. [T7 
reports an excess in P-vortex plaquettes of about 11%, at /? = 2.4, as compared to maximal 
center gauge. The greater number of P-vortex plaquettes found in Laplacian center gauge 
does not affect the asympotic string tension of projected configurations, and can probably 
be attributed to one (or both) of two sources: (i) a roughening of the P-vortex surface; 
(ii) some extraneous "small" P- vortices, unconnected to the single large vortex responsible 
for confinement, which is found [^ to percolate through the entire lattice. Either of these 
effects could modify the projected Creutz ratios at short distances, while preserving the 
asymptotic string tension. This may explain the slightly delayed approach of projected 
Creutz ratios to their asymptotic value, relative to previous results in maximal center 
gauge. 



Alexandrou et al. in ref. |T^ also suggest an alternative method for identifying vortices 
via Laplacian center gauge, which is based on locating certain gauge-fixing ambiguities, 
rather than center projection. We have not yet investigated this method; it would be 
interesting to know if our argument for the vortex-finding property in adjoint gauges can 
also be applied in this alternative framework. 

4 Conclusions 

Gauge-fixing has had a bad reputation, in connection with studies of the confinement 
mechanism. One is justifiably suspicious of any calculation which depends in an essential 
way on some special gauge choice, particularly if the physical motivation of that special 
gauge choice is unclear. In this article we hope to have dispelled some of that suspicion, at 
least in connection with maximal center and related gauges. Center vortices are created by 
discontinuous gauge transformations, which of course make no reference to any particular 
gauge condition. We have argued above that in maximal center gauge — and in an infinite 
class of other adjoint gauges — such discontinuous transformations are squeezed to the 
identity everywhere except on Dirac volumes, whose locations (together with those of the 
associated vortices) are then revealed upon center projection. This is the "vortex-finding 
property" which motivates the use of adjoint gauges, and it explains how adjoint gauge- 
fixing, combined with center projection, can extract the vortex content of thermalized 
lattice configurations. 

While we think it likely that every adjoint gauge, as defined in section 2, has the VF- 
property, this is certainly not true of every gauge-fixing procedure. The argument for the 
VF-property assumed a complete and unique adjoint gauge-fixing, but unfortunately the 
standard over-relaxation and simulated annealing methods are plagued by Gribov copies. 
This is a large loophole in the argument for vortex-finding, and in certain cases, e.g. the 
Kovacs-Tomboulis procedure, the Gribov copy problem is apparently severe enough to 
destroy the VF-property. 
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It is not clear, at present, why the Gribov problem destroys the VF-property only 
in some gauge-fixing procedures, but not in others. What does seem clear, however, is 
that the vortex-finding property, and the center dominance of projected configurations, go 
hand-in-hand. Maximal center gauge, asymmetric adjoint gauge, modulus Landau gauge 
with maximal center preconditioning, and Laplacian center gauge all have the VF-property, 
and all exhibit center dominance. Conversely, in the (i) Kovacs-Tomboulis procedure, (ii) 
modulus Landau gauge without preconditioning, and (iii) adjoint Coulomb gauge, the VF- 
property is lost, and in none of these cases is there center dominance in the projected 
configurations. 

We conclude with a sort of tautology: To find center vortices, one must use a procedure 
with the vortex-finding property. If the gauge-fixing -|- projection procedure doesn't have the 
VF-property, or if that property is destroyed by some modification (e.g. by Landau gauge 
preconditioning), then center vortices are not correctly identified on thermalized lattices, 
and center dominance in the projected configuration is lost. This fact does not call into 
question the physical relevance of P-vortices found by our usual procedure (which has the 
vortex- finding property); that relevance is well established by the strong correlation that 
exists between these objects and gauge-invariant observables. Ideally, potential problems 
due to Gribov copies could be avoided altogether by a gauge-fixing procedure which fixes to 



a unique adjoint link configuration. Laplacian center gauge fl^] appears to be an example 
of just such a procedure. 
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